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Abstract 

Results of recent searches for new physics in CP violation in charm decays from the BABAR ex- 
periment are presented. These results include a measurement of D" - D° mixing and searches for 
CP violation in two-body D° decays, a search for CP violation in the charm decays D ± — > K®K ± 
and Df — > K^K^, K^tt^, and a search for direct CP violation in the singly-Cabibbo suppressed 
—7- K + K^ir^ decays. These studies are based on the final dataset collected by BABAR at the 
PEP-II B factory at SLAC in the period 1999-2008. No evidence of CP violation is found in these 
charm decays. The measured mixing parameter yep = [0.72 ± 0.18(stat) ± 0. 12 (syst)]% excludes 
the no- mixing null hypothesis with a significance of 3.3cr. 
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1 Introduction 



The Cabibbo-Kobayashi-Maskawa (CKM) paradigm of CP violation (CPV) in the Standard Model 
(SM) has been tested by the BABAR and Belle experiments with high precision in many overcon- 
strained measurements. Nevertheless SM leaves many unanswered questions. CPV is one of the 
three Sakharov necessary conditions to generate the asymmetry between matter-antimatter (baryo- 
genesis) observed in the Universe. The measured CKM weak phase is unable to provide enough 
CPV to explain the observed baryon asymmetry. New CPV sources are needed from New Physics 
(NP) beyond the SM. 

At the B-factories important areas of search for NP are processes which are expected at low level 
in SM and which could be enhanced by NP. In these NP searches at low energies, charm physics 
plays currently an important and increasing role. The BABAR, Belle and CDF measurements of 
flavor mixing in the neutral D meson system pQ show evidence of D° - D° mixing at 1 % level. These 
results are in agreement with SM predictions [21 El H] and sets constraints on possible contributions 
from many NP models [5]. 

Recently the LHCb collaboration has reported a first evidence of CP violation in D° decays to 
K + K~ and 7r + 7r - [6]. This evidence has been confirmed by the CDF Collaboration [7j. Given the 
SM expectation that CPV in charm sector should be at the level of 1CP 3 (or lower) [2111], this CPV 
evidence was rather unexpected at the present experimental sensitivity. Marginally compatible 
with the SM expectations, this CPV may be a manifestation of NP or of significant enhancements 
of penguin diagrams in charm decays [H [101 [Til E2] • 

In this talk I present recent BABAR measurements concerning mixing and CP violation in charm 
sector. All new results presented here are preliminary. 

2 D°- D° Mixing and CP Violation in Two-Body D° decays 

Mixing in charm sector is unique because it involves virtual down type quarks. It arises from both 
short range and long range contributions. The short range contribution is expected to be very small 
because of CKM and GIM suppressions. The dominant long range contribution is non perturbative 
and hard to evaluate. This implies large theoretical uncertainties in the SM calculations of the 
mixing parameters x and y [SIEGE HI]- SM expectations values for x and y are < 10~ 3 but higher 
values are predicted in some NP models [151 US] ■ 

In a recent BABAR analysis [17] charm mixing and CP violation are measured using the ratio 
of lifetimes obtained in the D° decays to the two-body final states K^ir^, K~K + , and 7r~7r + . The 
analysis is based on an integrated luminosity of 468 fb^ 1 collected by the BABAR detector [IB] . 
Five different signal channels are considered [19]: three flavor tagged channels D* + — > D°Trf 
with D° —7- K + K~, D° — > it + it~, and D° — > K~ir + , K + ir~ and two flavor untagged channels 
D° — > K + K~ and D° — > K~ir + , K + tt~, where 7r+ is a slow pion track used in the tagging 
algorithm. 

The experimental observables yep |2U| , sensitive to mixing, and Al", sensitive to CPV , are 
measured. These observables are defined as: 

r + + f + r+ - f + 

y C p = f 1 and AY = — — — , (1) 

where r + (f + ) is the average width of the D° (D°) when reconstructed in the CP-even eigenstates 
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(K + K , 7T + 7r ). r is the average D° width describing the decays to the CP-mixed final states 

i^TT* PI]. 

Neglecting contribution of direct CP violation estimated at a level below our sensitivity [22] 
and taking into account that the CP violating weak phase (j) in SM to a good approximation does 
not depend on the final states [23] , the observables yc p and Ay in terms of the mixing parameters 
x and y can be written as: 

A i ^ M ■ A 
Ucp = ycos</>H — — xsm<p 

AY = —xsiiKp+'^^-ycoscj), (2) 

where Am = (\q/p\ 2 — \p/q\ 2 )/(\q/p\ 2 + \p/q\ 2 ) measures the CP asimmetry in mixing. The complex 
parameters p and q relate the mass eigenstates of neutral mesons , \D% i >, to the flavor eigenstates, 
\D° > and \D° >, through the relation |Z>i ;2 >= p\D° > ±q\D° >. 

A simultaneous extended unbinned Maximum Likelihood (ML) fit to the two-dimensional dis- 
tribution of the proper time and proper time error in tagged and untagged modes is performed: the 
average D° lifetime r is extracted from K^n^ final states and the effective lifetime r + (f + ) is ex- 
tracted from D° (D°) decays to the final states K~ K + and ir~ir + . Main sources of background are 
misreconstructed charm events and the combinatorial background candidates consisting of random 
tracks. Using the reciprocals of the three measured lifetimes in Eq. [T]we obtain: 

ycP = (0.72 ±0.18 ±0.12)% and AY = (0.09 ± 0.26 ± 0.09)% , 

where the first uncertainty is statistical and the second systematic. Projections of lifetime fit are 
shown in Fig. [T] 

These results exclude no- mixing hypothesis at 3.3<r significance and show no evidence of CPV . 
The yep value is consistent with the mixing parameter y measured in the decays D° — > K®h + h~ 
(where h = K,ir) [23] as expected in absence of CPV. This yep measurement is the most precise 
single measurement up to date. These results are in agreement with SM predictions. 



3 CP Violation in D ± -4 and Df -> K"^, 

The channels D ± -»• can proceed through Cabibbo-Favored (CF) and Doubly-Cabibbo- 

Suppressed (DCS) transitions. The CF transition is largely dominating and the SM expectation 
for direct CP is negligible. The channels Djr — > K^K^ , K^ir^ can proceed through two Singly- 
Cabibbo-Suppressed (SCS) transitions, both of comparable amplitudes. The relative phase between 
the two decay amplitudes can generate interference effects and induce direct CPV . 

In these channels with a K° (or K ) in the final state a time-integrated CPV of ~ (±0.332 ± 
0.006) % is induced by the K°K° mixing [25J. The sign of this asymmetry is positive (negative) in 
presence in the final state of a K° (K°). The exact value of this CPV asymmetry contribution de- 
pends on the requirements on the reconstructed K® — > 7r + 7r~ decays and the decay kinematics |26j. 

Previous results of searches for direct CP violation in these decay modes by CLEO-c [27] and 
Belle [28] Collaborations are all in agreement with SM expectations. 

Direct CP asymmetry in these charm decay modes has been recently searched for by BABAR [2U] 
using a dataset of 469/6 -1 . The following direct CP- violating parameter Acp is measured for each 
decay channel: 
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Figure 1: Proper time fit projections with the fit results overlaid. The combinatorial background 
distribution (Comb) is stacked on the top of the misreconstructed-charm background distribution 
(Charm). Under each plot are shown the normalized Poisson pulls; "unt" refers to the untagged 
dataset. The gray band is the PDG D° lifetime ±<t|25|. 



_ r(D+ -> K°(n+,K+)) - T(D-(s) -> K°(ir+,K+)) 
CP ~ r(D+ s) ^K0(7r+,K+)) + r(D+ s) ^K°(ir+,K+)) ' 

where T is the partial width of the decay channel. 

The measured CP asymmetry A an be written as A = Acp + Am + A e , where Acp is the direct 
CP asymmetry contribution, Afb is a forward/backward asymmetry contribution in cc production 
from 7 — Z° interference and higher order QED processes, and A e is the asymmetry contribution 
induced by the detector in tracking, particle identification, and in material interactions. Afb 
asymmetry is an odd function of the cosine of the polar angle of the meson momentum in the 
e + e~ center of mass (CM) system, cosO^. Acp and Afb are both measured while data have been 
corrected for A e with a control sample. The data-driven method used to correct for A e is described 
in Ref.|3DJ. 

A simultaneous binned ML fit to the and D7 s s invariant mass distributions is performed 
in 10 equally spaced bins of cos6* D with bin at [-1.0, -0.8]. Since Acp is independent of the 
kinematic variable cos^, the asymmetry -4(+| cos 9* D \ ) measured in a positive cos9* D bin and 
the asymmetry A{—\cos9* D \ ) measured in its symmetric (negative) counterpart cos6* D bin give 
the same contribution to Acp. On the other hand since Afb is an odd function of cos6* D , the 
contribution to Afb from symmetric cosO^ bins have the same magnitude and opposte sign. So 
Acp and Afb as a function of cos 6^ can be written in the form: 
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A FB (\cos* e D \) 
A C p{\ cos* 6 D \) 



A{+\cos6* D \) - A{-\cos6* D \ 
2 

A{+\cosO* D \) + A(-\cosQ* D \ 

2~ 



The values of Acp and Afb asymmetries are shown in Fig. [2j 
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Figure 2: .4 CP (top) and .Aps (bottom) asymmetries for D ± ->■ -fi^-K^ (left), -> 
(center), and .D^ — )• -f^gTr 1 * 1 (right) as a function of | cos* #d| in the data sample. The solid line rep- 
resents the central value of Acp and the gray region is the ±cr interval, both from a \ 2 minimization 
assuming no dependence of Acp on | cos* 9d\- 



For each decay mode Table [T] shows the Acp value from the fit, the bias corrected Acp value 
and in the last raw the Acp value after subtracting the expcted Acp contribution due to K°-K° 
mixing. These results are consistent with zero and with the SM predictions within 1 a. 



4 CP Violation in the Decays D 
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Searches for direct CP violation in the SCS decays — > K+K^tt^ have been recently performed 
by BABAR, using a dataset of 476/6~ 1 |31j. This sample contains enough 3-body SCS decays 
to probe CP at the level of SM predictions. The decay D + — > K + K~tt + |19| is dominated by 
quasi-two body decays with resonant intermediate states, giving possibility to study direct CPV 
in a particular resonance or in different regions of the Dalitz plot (DP). In previous analyses of 
these 3-body decay modes performed by CLEO-c [32] and LHCb [33] Collaborations no evidence 
for CPV has been found in agreement with SM prediction. 

Signal reconstruction efficiency is determined with a sample of Monte Carlo (MC) simulated 
events from the distribution of reconstructed events as a function of the CM polar angle of the D 
meson (cosOcm) and of the m 2 (K~Tr + ) vs m 2 (K + K~) DP. The ratio of efficiency-corrected signal 
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Table 1: Aqp measurements. Uncertainties, where reported, are statistical the first and systematic 
the second (BABAR Preliminary). 





D ± ->■ K"K ± 


Df KgK* 


Df -> 


Acp value from the fit 


(+0.16 ±0.36)% 


(0.00 + 0.23)% 


(+0.6 + 2.0)% 


Bias Corrections for: 
Toy MC experiments 
PID selectors 
Ks — K® interference 


+0.013% 
-0.05% 
+0.015% 


-0.01% 
-0.05% 
+0.014% 


-0.05% 
-0.008% 


Acp final value 


(+0.13 + 0.36 + 0.25)% 


(-0.05 + 0.23 + 0.24)% 


(+0.6 + 2.0 + 0.3)% 


Acp contribution 
from K° — K° mixing 


(-0.332 + 0.006)% 


(-0.332 + 0.006)% 


(+0.332 ± 0.006)% 


Acp final value (charm only) 


(+0.46 + 0.36 + 0.25)% 


(+0.28 + 0.23 + 0.24)% 


(+0.3 + 2.0 + 0.3)% 



yields, R = ^ T D+ ^ D+ = 1.020 ± 0.006 is used to allow for asymmetries in the MC event production 
due to physics or detector induced effects. 

Time-integrated CP asymmetry (charge asymmetry) is defined in a given bin as: 

A ^ N D+/^- N D-/^- (5) 

N D+ /e D+ +N D -/e D - 

Selection efficiencies are corrected to account for differences between data and MC simulated 
events in the reconstruction asymmetry of charged pion tracks and in the production model of 
charm mesons. The charge asymmetry A contains contributions from both the forward/backward 
asymmetry Afb an d the direct CP asymmetry Acp- To remove the contribution Afb, the charge 
asymmetry is averaged over four symmetric bins in cos 6 cm- The averaged values of Acp in the 
four bins are shown Fig. [5J The central value Acp = (0.35 ± 0.30 ± 0.15) % is obtained with a \ 2 
minimization. The probability that the asymmetries are null in all the four bins is 21%. 

A model-independent technique to search for CP violation in DP is to compare CP asymmetries 
in different regions of the DP. The results of Acp asymmetry measured in four regions of DP are 
given in Table [2j Measured CP asymmetries are consistent with zero. 

We also measure the normalized residuals A of efficiency-corrected and background subtracted 
DP for D + and D~ for equally populated bins. A is defined as 

^ n{D+)-Rn{D-) 
v / a 2 (D+)+R 2 a 2 (D-) ' 

where n is the yield in a bin in the DP and a its uncertainty. 

A distribution is fitted to a Gaussian function. For 100 bins we obtain a Gaussian residual 
mean of 0.08 ± 0.15 and a width of 1.11 ± 0.15. The probability that the two DP's are consistent 
with no CP asymmetry is 72%. 

Angular moments of the cosine of the helicity angle 9h of the D decay products reflect the 
spin and mass of intermediate resonant and non resonant states [M]. The helicity angle 6h in 
the decay D — ¥ (r — > AB)C is defined as the angle between the momenta of B and parent D in 
the AB rest frame. We can search for CP in the DP in a model-independent way by comparing 
the angular moments between D + and D~ [35]. Angular moments of order 1 are defined as the 
efficiency-corrected and background-subtracted two-body invariant mass distributions (m(K + — , 
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Figure 3: Charge asymmetry Acp as a function of | cos 9cm I m data. The solid line represents 
the central value of Acp and the dashed lines represent the ±cr interval, both determined from a 
X 2 minimization assuming no dependence on | cos 0cm I • 



Table 2: Yields, efficiencies,, and CP asymmetry in four different regions of the DP. First uncer- 
tainty in the CP asymmetry is statistical, the second is systematic (BABAR Preliminary). 



Dalitz plot region 


N(D + ) 


e(D+)[%] 


N(D-) 




A C p[%] 


Below ir(892) u 


1882 ± 70 


7.00 


1859 ± 90 


6.97 


-0.65 ± 1.64 ± 1.73 


K* (892)° 


36770 ± 251 


7.53 


36262 ± 257 


7.53 


-0.28± 0.37 ± 0.21 


0(1020) 


48856 ± 289 


8.57 


48009 ± 289 


8.54 


-0.26 ± 0.32 ± 0.45 


Above K* (892)° and 0(1020) 


25616 ± 244 


8.01 


24560 ± 242 


8.00 


1.05 ± 0.45 ± 0.31 



m(K vr + )) weighted by spherical harmonic moments = Y^cosOh). Weights in two-body 
invariant mass intervals are defined as: 



W f) = V J L ; (7) 

v. ti s> 

where i is bin index, and j, k event indices. S and B refer to signal and background, and < ej > 
is the average efficiency in bin i. 

Normalized moment residuals Xi for D + and D~ are calculated for 1 from to 7: 

(wP(d+)-rw«\d-)) 

/<t<P\d+) + B*v<P\d-) 
The x 2 is calculated over all the mass bins in K + K~ and K~ir + moments with 

i L U 
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Table 3: CPV parameters from the DP fit. First uncertainties are statistical, the second systematic 
(BABAR Preliminary). 



Resonance 


r (%) 


A</> (°) 


is* /con\0 
A (OilZ) 


U. (bLJVrjD) 




if* 1 1 A QH"\0 

A (1430) 


O -i pi + 5 . 6 5 i a a o 
- 9 - 4U -S.36 ± 4 - 42 


f? i i +3.29 i -t qo 
— D.ll_ 3 24 ± l.oy 




(1 QCl+O- 82 _i_ n en 


7 /1Q+ 3 - 55 ±9 Qt; 


NR 


-14.30±^;|? ±5.98 


-2.56±6;°7±8.91 


k(800) 


2.00±|;™±1.85 


2.10tl\f 5 ± 1-01 


a (1450)° 


5.07te 54 ± 9.39 


4.00^3-96 ± 3.83 




Ax 


Ay 


/o(980) 


-0.199+o i°o ±0.084 


-0.231+o;i°° ± 0.079 


/o(1370) 


0.019l°;o48 ± 0.022 


-0.0045l° o39 ± 0.016 



where p\ 1 2 is the correlation coefficient between X\ 1 and X\ . 

With a number of degrees of freedom NDF equal to 287 the x 2 /NDF in the K + K~ and K~ir + 
moments is 1.10 and 1.09, consistent with no CPV at 11% and 13 %, respectively. 

BABAR also searched for CPV in a model-dependent DP analysis of the D + — > K + K~ir + de- 
cay |19j . The DP amplitude A in the isobar model is written as a set of two-body intermediate 
states r : A = ^ r M. r e l ^ r F r , where M. r and 4> r are real and F r = F r (m(K + K~),m(K~ir + )) are 
dynamical functions describing the intermediate states. In case of amplitudes with small contri- 
butions the complex coefficient has been parameterized in a Cartesian form: x r = M r cos(fc r and 
y r = M. r sin (ft r . The K*(892)° has been chosen as the reference amplitude. Assuming no CPV 
the relative fractions of resonances and a constant non resonant amplitude over the entire DP 
contributing to the decay have been determined with an unbinned ML fit. 

To allow for possible CPV in the decay, the resonances of the D + (D~) decays contributing 
with a fit fraction of at least 1% have been parameterized with different amplitudes and phases 
in their decay amplitudes. A simultaneous fit to the D + and D~ samples have been performed, 
parameterizing each resonance with four parameters, M r , (fi r , rep, and A<pcp- The CPV parameters 
are rep = |^ r | 2 ~|^ r | 2 and A(f>cp = <j) r — 4> r ■ The Cartesian form of the CP violating parameters 
are Ax r and Ay r with x r (D^) = x r ± Ax r /2 and y r (D^) = y r ± Ay r /2. Fit results are shown in 
Table [3 

All CPV parameters from DP fit are consistent with zero and with SM expectations. 

5 Conclusions 

I have presented recent improved BABAR measurement of mixing and a search of CP violation 
in two-body D° decays, a search for direct CP violation in D^ — > K®K^ and D^r — > K^K^ 1 , 
Kgir^, and searches for CP violation in the decays D^ — > K + K~ir ziz both using model-independent 
and model-dependent analysis techniques. All results in these analyses are well described within 
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the SM and no effect related to NP has been found. The measured mixing parameter yep = 
[0.72 ± 0.18(stat) ±0.12(syst)]% excludes the no- mixing null hypothesis with a significance of 3.3o". 
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